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Summary 

Equations for the semianalytic representation of 
a class of surfaces that vary smoothly in cross- 
sectional shape are presented. Some methods for fit- 
ting together and superimposing such surfaces are 
described. A brief discussion is also included of the 
application of the theory in various contexts such 
as computerized lofting of aerodynamic surfaces and 
grid generation. 

Introduction 

A very general class of surfaces can be approx- 
imated for computational purposes by a set of fiat 
panels or of curved patches (ref. 1). Aside from the 
storage and accounting problems of defining and re- 
taining the information for each panel, such repre- 
sentations present difficulties for applications that 
require continuity of derivatives at the panel edges. 
The more stringent the smoothness requirement at 
the edge, the higher the order of the polynomial that 
is required for the modeling. This type of polynomial 
modeling leads to a surface “waviness” that may be 
difficult to detect visually but manifests itself in an 
oscillation of the derivatives. 

This problem does not exist for simple config- 
urations that can be synthesized from basic geo- 
metric shapes for which exact analytic descriptions 
are known. However, for aerodynamic configura- 
tions, this limitation to basic shapes is too restric- 
tive for meaningful design. The present paper ex- 
plores a somewhat different approach — one that re- 
tains many of the advantages of fully analytic rep- 
resentation while enlarging the class of shapes rep- 
resentable beyond simple geometric objects to in- 
clude aerodynamic-type surfaces. A brief descrip- 
tion is included of some advantages of the method 
relative to several problems: computerized lofting 
of aerodynamic shapes (both external and internal) 
and gridding for flow-field calculations. The proce- 
dure is a generalized version of a method described in 
reference 2. 


Symbols 


A 

surface area 

c 

shape transition function 

E,F,G 

metric coefficients 

i,j,k 

orthonormal base vectors 

m 

exponent parameter (eq. (5)) 

N 

unit surface normal vector 

r 

surface vector 


R, 9 polar coordinates 

t independent variable used to define 

surface cross-sectional shape 

t/ parameter defined by equation (29) 

x,y,z Cartesian coordinates 

x m = 1 - x 3 

■q independent variable used for two- 

dimensional grid generation 

A size transition function 

fi homotopy function of £ 

£ independent variable used for three- 

dimensional grid generation 

Subscripts: 

1,2,3 various cross-sectional shapes used in 
defining transition surface 

a,b,c quantities associated with different 

surfaces 

/ fuselage surface 

i inner boundary 

l limiting 

n normalized 

o outer boundary 

s canopy surface 

t, x derivative with respect to t or x, 

respectively 

An asterisk (*) with a symbol denotes a c- function 
used in grid generation. A tilde (~) over a symbol 
denotes a A-function used in grid generation. Primes 
denote a derivative with respect to the argument. 

Analysis 

Surface Equations 

The analysis is described as it applies to the ba- 
sic lofting problem in its simplest form: a surface 
component is to be designed so that it varies gradu- 
ally and smoothly from a given initial cross-sectional 
shape to a different specified terminal or base shape. 
Such a component might represent, for example, a 
forebody, a wing, or a duct. The direction of varia- 
tion is taken to be the x-direction, with the two end 
shapes specified in y-z planes (fig. 1(a)). 

Neglecting for the moment the variation in cross- 
sectional size, consider first only the shape variation. 


(8b) 


Suppose the initial shape can be represented in para- 
metric form with parameter t: 


3/i = Vi(t) 1 
zi = zi (0 / 


{t a <t< t b ) (1) 


and the terminal cross section is also representable 
in parametric form with the same t domain (t a ,t b ): 

2/2 = 1 / 2(0 (2a) 

z 2 = z 2 {t) (2b) 

(See fig. 1(b).) 

For example, if the cross sections are appropri- 
ately represented in polar form, then 


y 1 = R^fljcos# 
zi = Ri(0) sin 0 


(0 < 0 < 2tt) 


with similar expressions for y 2 and 2 2 . The x-variable 
is now normalized by 


X — Xj 

2:2 - Zl 


(3) 


and a parameter c n (x n ) is defined such that c n varies 
from 0 to 1 as x n varies continuously and monotoni- 
cally from 0 to 1 (i.e., a homotopy parameter). Then 
define c by 

c(x) = c n [x n (x)] (4) 

Thus, c(xi) = 0 and c(x 2 ) = 1. The function c(x n ) 
is usually specified to be smooth and monotonic. A 
typical expression for c n might be 


c„(x n ) = x™ (5) 

with m as an adjustable design parameter. A tran- 
sition surface representing a smooth blending of the 
two end shapes is defined by the functions 

y{x,t) = [1 - c{x)\yi (f) + c(x) y 2 (t) (6a) 

z(x, t) = [1 - c(x)] 2 !(£) + c(x) z 2 {t) (6b) 

Finally, the size variation of the surface is pre- 
scribed by introducing an independent scaling func- 
tion A(x) which is specified to be smooth and non- 
negative but not necessarily monotonic. Thus, the 
lofted surface is defined by the equations 

y{x, t) = A(x){[l - c(x)]y x (t) + c(x) y 2 {t)} (7a) 

z{x,t) = A(x){[l - c(x)\z x {t) + c(x) z 2 {t)} (7b) 

A somewhat different type of transition surface is 
defined by 

2 /(x, t) — Ay| lc) 2 / 2 (8a) 


z{x,t) = Xz[ 1 c) y 2 

where c and A are functions of x, and y x , y 2 , z x , and 
z 2 are all functions of t. Other, more complicated 
forms of transition surface equations may also be 
defined. The vector form of the surface equation is 

r(x,t) = xi + y(x,t)j + z(x,t)k (9) 


Computation of Surface Parameters 

Equation (9) together with equation (7) or (8), 
defines the transition surface in terms of four func- 
tions, each defined as a function of a single inde- 
pendent variable. The two end shapes ( 2 / 1 , 21 ) and 
{V2,z 2 ) depend only on t, whereas the shape varia- 
tion function c and the size variation function A de- 
pend only on x. Therefore, the derivatives of the 
surface vector r(x, t) can be computed in terms of 
the derivatives of these individual functions. Thus, 

r x (x, t) = i + [(A x /A) y(x, t) + A c x (y 2 - y x )\j 

+ [(A x /A) z(x,t) + \c x (z 2 - 2 i)]fc (10) 

and 

r t (x, t) = Of + A[(l - c)y lt + cy 2t ]j 

+ A[(l - c) 2 i t + cz 2t ]k (11) 

(See fig. 1(c).) In equations (10) and (11), c and A are 
functions of x, while 2 / 1 , 2 / 2 ) z i , and z 2 are functions 
of t. With these equations, one can directly compute 
the surface parameters that depend only on the first 
derivatives, such as the metric coefficients 


E = r x ■ r x 

(12a) 

F = r x r t 

(12b) 

II 

n 

(12c) 

the surface normal 


n - T 1 x rt _ 

y/EG - F 2 

(13) 

and the local area element 


dA = VEG - F 2 dx dt 

(14) 


By differentiating equations (10) and (11) again, the 
second derivatives r xx , r xt , and r t t can be obtained. 
From these, the local curvature coefficients can be 
calculated as well as surface parameters such as the 
principal curvature, mean and total curvature, and 
curvature directions according to the basic differen- 
tial geometry formulas. (See, for example, ref. 3.) 
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Computed Example 

An example of a transition surface is shown in 
figure 2. This surface, which might represent the 
upper surface of a blended wing-body combination, 
was generated by equations (7) and (9) with the 
following choice of variables. For the initial cross 
section at aq = 0, 

yi(t)=t \ 

v/T^P (-1<<<1) 
= 3 — J 

For the final cross section, at x 2 — 5.0, 

2/2(0 =t (-1 < t < 1) 

z 2 (t) = \J (0.317) 2 — f 2 (10 < 0.317) 

z 2 (t) = -O.2(|0 - 0-317) (1 > |0 > 0.317) 

The ^-coordinate was nondimensionalized as follows: 


x 



A linear variation of the shape transition was pre- 
scribed as follows: 

c n{ x n) = x n 

but, for this particular example, a more general 
version of the size variation function 

A(x n ) = A 2 (1.6x n - 5.316x2 
+ 9.133x^ _ 4.417x4) 

was required to provide the smooth transition surface 
shown in figure 2(b). 


can then be used in the transition surface equations 
(eqs. (7) to (11)), with the other end shape specified 
either analytically or numerically. 

Similarly, the shape- and size- variation functions, 
c(x) and A(x), could be specified numerically. How- 
ever, since these variations are normally smooth and 
gradual for aerodynamic shapes, analytic expressions 
are normally satisfactory for their representation. 

It should be emphasized that, if one or more of 
the basic geometry functions is defined numerically 
but smoothed so that the required derivatives are 
meaningful, then the surface parameters can still 
be computed from the transition surface formulas 
(eqs. (7) to (9)). This formulation, which represents 
an analytic combination of numerical and analytic 
functions, is termed “semianalytic.” 

Specifying an Intermediate Cross-Sectional 
Shape 

A designer may need to define a surface which 
is to vary gradually in cross-sectional shape but in 
such a way that intermediate cross sections cannot 
reasonably be described as a blending of the two end 
shapes. Two approaches to this problem are dis- 
cussed. The first approach is to define a single transi- 
tion surface that includes the specified intermediate 
cross-sectional shape. Denote this cross section at 
intermediate axial station x 3 by y 3 (t), z 3 (t) with t 
having the same domain as for the two end shapes. 
For the sake of simplicity in notation, x is normalized 
so that 

Xi = 0 

x 2 = 1 
0 < x 3 < 1 


Use of Digitized Input 

In the event that one of the specified cross- 
sectional end shapes is difficult to describe analyt- 
ically, it is necessary to represent it numerically (for 
example, by means of digitized data). In this case, 
the data must be smoothed in order to avoid the ir- 
regularity and oscillation of the derivatives that are 
characteristic of unsmoothed data. This smoothing 
is accomplished through the use of 6-splines. (De- 
tails on 6-spline theory and techniques are given in 
ref. 4.) The numerical cross-sectional data are dis- 
played graphically along with the first two numerical 
derivatives (fig. 3(a)). The various 6-spline param- 
eters are adjusted interactively until a shape is ob- 
tained that represents a satisfactory smooth approxi- 
mation to the original input shape and its derivatives 
are free of spurious oscillation (fig. 3(b)). These nu- 
merical functions describing the cross-sectional shape 


Define 


x m = 1 — 2:3 


Then the simplest surface that satisfies the pre- 
scribed conditions is given by 


y{x,t) = A(x)| 


, 1 + z 3 

1 x 


2/i (0 - 


x 2 ] 

H 

x 3] 

+ - — ~ x ysit) 

x 3 x m 


x 3 

x 3 


' x 2/2 (0 


(15) 


with a similar expression for z(x,t). 

This formula satisfies the required conditions, but 
it contains no free parameters (other than A(x)) to 
permit further flexibility in design. Some freedom in 
design could be introduced by the addition of higher 
degree terms in x, but at the expense of greater 
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complexity and with the risk of creating a waviness 
in the surface. 

The second approach for introducing a speci- 
fied cross-sectional shape between the two given end 
shapes is to treat the surface as two transition sur- 
faces joined end to end at x 3 . Define the new nor- 
malized axial variables as 


If the second derivative were required to be con- 
tinuous at x 3 , then in addition to equation (22), the 
condition 

c"(l) = c"(0) = 0 (24) 

would be enforced. 

If a fourth cross section is specified at x 4 (where 
£3 < x 4 < x 2 ) so that a third transition surface (with 
transition function c c ) between x 3 and x 4 is required, 
then it must satisfy the conditions 


x — x 3 

x b = 

x 2 - x 3 

Since multiplication by a smooth scale function A(x) 
will not affect the differentiability of the resulting 
surface, the two unsealed transition surfaces will be 
considered: 



ya{x a ,t) = (1 - c a )y l + c a y 3 (17a) 

Z a{x a <, f) (1 ^a)^l + CaZ 3 (17b) 

y b {x b ,t) = (1 - c b )y 3 + c b y 2 (18a) 

z b {x b ,t) = (1 - c b )z 3 +c b z 2 (18b) 

At the juncture location x 3 , 


c a = 1, c b - 0 
x a — 1, x b - 0 


Va = V 6 


2/a, = y b . 


2/a* = (ys - y i) 


dCg dXg 

dxg dx 


yz - y i 

X 3 — Xj 


(19) 

( 20 ) 

(21a) 


y bl = 


V2-V3 , 
c 6 

X2 ~ X 3 


(21b) 


with similar expressons for the 2 -coordinate. In 
equations (17), (18), and (19), c a and c b are functions 
of x, whereas y\, y 2 , y 3 , X\, x 2 , and x 3 are functions 
of t. 

From equations (21), it is apparent that the first 
derivative with respect to x is continuous at x 3 only 
if Cg(xg) and c b (x b ) are defined so that 


4 ( 1 ) = 4 ( 0 ) = 0 ( 22 ) 


except in the unusual circumstance that the two 
functions y 2 (f) and yi(t) are identically equal, in 
which case the requirement is that, at x 3 , 


= r4(°) < 23 > 

x 3 ~ x 3 

For this equation to be satisfied, one of the shape 
transition functions would fail to meet the mono- 
tonicity condition. 


4(0) = 4(1) = 0 (25) 

in order to ensure continuity of the first derivative at 
both ends of the interval. 

In theory, this procedure could be continued in- 
definitely with the insertion of an arbitrary number 
of specified intermediate cross sections. However, 
the requirement that the shapes be “lined-up” at 
the juncture according to equation (25) would tend 
to lead rapidly to the undesirable waviness of the 
surface. 

Combining Surfaces 

After each individual component has been rep- 
resented in analytic or semianalytic form, there re- 
mains the problem of synthesizing these components 
into a configuration. Discussion of all the various 
approaches to this problem is beyond the scope of 
the present paper, but one technique that has been 
used effectively is simple superposition. Its applica- 
tion in the present context is illustrated with a partic- 
ular example — that of superimposing a canopy onto 
a fuselage (fig. 4(a)). Let the equation for the upper 
surface of the fuselage be 

*•/(*,«/) = xi + y/(x, t/)j 

+ Zf(x,tj)k (0 < x < x 2 ) (26a) 

and let the equation for the canopy surface be 
r »(x,f s ) = xi + y 9 {x,t a )j 

-I- z 3 (x, t 9 )k (0 < xi < x < x 3 < x 2 ) 

(26b) 

For purposes of illustration, assume that the range 
of t is (— 1, +1) for both 17 and r s . Then t = 
— 1 and t = +1 represent the lateral extremes of 
the components, which are symmetric about t = 
0. Let Xj denote the axial station at which the 
canopy surface is to be initiated. Then, for any 
x > Xj, y s {x — Xj,l) is calculated. Then, for this 
value of y, the corresponding value of t for the 
fuselage surface is calculated (fig. 4(a)). Suppose, 
for example, that 

y/(x, tf) = Xf(x)if (27a) 
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y a {x-xi,t s ) = X s (x-x 1 )t a (27b) 

y s (x-xi, 1) = A s (x-xi) (28) 

Then, for the fuselage surface, the corresponding 
value of tf is found by substituting relation (28) into 
relation (27a) 

_ A a (x -i!) 

i X f( x ) 

Such a boundary value, £*, can be determined for 
each value of x > X\. If an analytic expression 
such as equation (29) cannot be obtained, ti must be 
computed numerically. Then the synthesized surface 
is defined by the relations 

r(x,tf) = Tf(x,tf) (x < Xi) (30a) 

r(x,t/) = Tf(x,tf) (x > in \tf\ > |t*(x)|) (30b) 

r(i, t s ) = x i + y s {x-x \,t s )j 
+ [z/(x,t f ) 

+ x s (x-xi,t s )]fc (x>xj; |t/| < |*z(x)|) (30c) 

Figure 4(b) shows a surface formed by a synthesis of 
several components in this manner to form a fuselage- 
like geometry. 

Discussion Regarding Applications 

Aerodynamic Design 

The transition surface equations (eqs. (7) to (9)) 
were originally developed as a means for rapid calcu- 
lation of surfaces generated by a lofting procedure. 
The two end cross sections are specified by the de- 
signer, and the intermediate cross sections are com- 
puted from equation (9) as x = Constant lines. These 
can be displayed graphically; and the shape and scale 
functions, c(x) and A(x), can be varied interactively 
to obtain the desired transition. The lofting lines are 
the t = Constant lines. Since the two end shapes 
are specified so that there is a natural one-to-one 
t-correspondence, this correspondence holds through- 
out the transition, and consequently there is no prob- 
lem with lofting lines developing waves or crossing 
each other. 

Grid generation 

A single transition surface such as that shown in 
figure 2(b) is especially easy to grid for purposes of 
flow calculations. The (x, f ) coordinate lines form a 
natural surface grid. The external field can be grid- 
ded by a simple extension of the homotopic procedure 
that is used to generate the surface itself. 

For calculation of supersonic flow over such a 
surface, a quasi-two-dimensional grid is sufficient. 



An outer boundary surface must be specified. If the 
configuration is contained within a duct (e.g., a wind- 
tunnel test section) , the duct wall serves this purpose, 
but otherwise a simple shape, such as a cylinder, 
is arbitrarily chosen. This outer boundary shape 
is also represented as a transition surface with the 
same distribution of the t and x variables as used to 
describe the inner surface. Then the grid is generated 
by defining, at each constant x-station, a homotopic 
variation from the inner surface t/i, z x to the outer 
boundary y oy z 0 . Denote this homotopy parameter 
by rj, and the shape transition function by c*. Thus, 

c*(rj) = 0 (r, = 0) (31a) 

c*(r?) = l (*7 = 1) (31b) 

This function is taken to be independent of x; it is 
the same at each x-station. 

Now denote the size variation function by A* so 
that < A* < A c . Since the size of the inner surface, 
and possibly the outer boundary also, varies with x, 
A* cannot be independent of x. It can, however, be 
expressed in terms of a normalized function A by 

A* = A t + A(A 0 -A t ) (32) 

so that 

A(*7) = 0 (rj = 0) 

A (v) = 1 (*7 = 1) 

and A is independent of x. By defining the grid 
at each x-station in terms of these functions, c*(r?) 
and A(r?), a natural correspondence is established 
between grid points at the various x-stations. If 
both inner and outer boundaries vary smoothly with 
x, the corresponding grid-point locations will vary 
smoothly. Figure 5 shows a grid generated in this 
manner for a cross section of a transition surface 
similar to that of figure 2(b). 

For some flow calculations, such a quasi-two- 
dimensional grid is not always adequate. If the 
grid structure is required to reflect the influence of 
the body shape upstream of the nose, a fully three- 
dimensional grid can be defined as follows. Take the 
average x-station for the surface as the origin 

so that the surface is now defined on the domain 

-a* < e < e.,i (-!<*< i) 

For the outer boundary shape, define an appropriate 
convex surface (e.g., ellipsoid) on the same £, t do- 
main with the same distribution of the £ and t vari- 
ables as for the inner surface. The outer boundary 
shape is obtained by multiplying this surface r(£, t) 
by a large scale factor A 0 . The grid homotopy pa- 
rameter is again denoted by r/; the shape transition 
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function, by c*(r?); and the size transition, by A* (77) 
(where 1 < A* < A 0 ). Also a homotopic function /z(£) 
is defined such that /z varies monotonically from 0 
to 1 as |f | varies from 0 to 1. Then the grid-point 
locations are given by 

= A*(?7) /z(f) (33a) 

y(Z,t,v) = A*(r?){[l - C*(r?)]y,(f,t) 

+ c*{r]) y a {r),t)} (33b) 

= A*(r?){[l - c*(Tj')]z i (f,<) 

+ c*(ri) z 0 (£,t)} (33c) 

As in the quasi-two-dimensional case, considerable 
flexibility is permitted in designing the grid by in- 
dependent control of the shape and size variation 
functions. 

A grid-type structure generated by equations (33) 
is illustrated in figure 6. Figure 6(a) shows a tapered, 
cambered wing. Figure 6(b) shows two of the f = 

Constant sections viewed from the wing root. Fig- 

ure 6(c) shows a t — Constant section viewed from 
the front of the wing. A blow-up of the circumscribed 
region in figure 6(c) is shown in figure 6(d) to display 
the shape of the f = Constant lines near the surface. 

Concluding Remarks 

Equations for the semianalytic representation of a 
class of surfaces that vary gradually in cross-sectional 


shape have been presented. It was demonstrated that 
basic surface geometry parameters can easily be calcu- 
lated with such a representation. Some methods for 
fitting together and superimposing the surfaces were 
described. Also included was a brief discussion of the 
application of the theory in various contexts such as 
computerized lofting of aerodynamic surfaces and grid 
generation. 

NASA Langley Research Center 
Hampton, VA 23665 
December 18, 1984 
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Station 2 
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(a) Perspective view of lofted surface. 



(b) Normalized end sections. 
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(a) Basic surface: tapered, cambered wing. 

Figure 6. Sample three-dimensional grid-type structure for transition surface 
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